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Abstract — In this study, the frictional losses in a pipe due to 
shear stress produced in a pipe due to viscosity of fluids has been 
discussed. Pipeline system is a vast network in any industry and 
minimizing losses through these pipes can contribute a lot. This 
paper contains flow analysis of different fluids in a pipe 
entrance. The Computational Fluid Dynamics (CFD) analysis 
for the flow of non-Newtonian fluid through a straight pipe is 
presented. Laminar non-Newtonian pseudoplastic power law 
model has been used for the simulation of non-Newtonian fluid 
flow through the pipeline, particularly the entrance area. The 
CFD analysis has been tested from previously published 
experimental results. 


Index Terms — CFD, Pipeline, Channel, Non-Newtonian 
fluid, Pressure drop. 
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I. INTRODUCTION 

Implementation of fluids in many industrial activities often 
requires their transport in pipeline networks having 
singularities, pumps and heat exchangers. A good design of 
these networks requires a precise knowledge of the laws 
governing the transfer of heat energy loss E (regular and 
singular). 

Fluids with nonlinear [3] behaviour conforming the power 
law as well as their changing flow parameters are of interest to 
many industrial sectors; they are nevertheless considerably 
less studied compared to Newtonian fluids. The consequences 
of shear thinning need to be identified to a clearer 
understanding of the structure of the flow. 

The assumption that the flow is fully developed in regions 
where it remains under strong influence of the inlet boundary 
conditions can seriously underestimate the design of flow 
systems and incorrectly assume specific velocity profile 
shapes leading to wrong conclusions in the interpretation of 
data [10]. 

Unsurprisingly, determining this length for a range of ducts, 
e.g., pipe, channel, and annuli, has been the subject of a great 
deal of attention over the past 100 years or so [9]. 

A detailed discussion of the inconsistencies and confusion in 
literature is provided by Durst et al. [9], who conducted a 
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detailed numerical study and proposed nonlinear correlations 
for pipes and channel for Newtonian fluids. 

In this study we are interested in the laminar flow of 
non-Newtonians pseudoplastics thermo-dependents fluids in 
a straight cylindrical pipe and between two plates particularly 
in the entrance area. 

The mathematical model is represented by the equations of 
continuity and conservation of momentum. The rheological 
behaviour is given by the power law of Ostwald of Wale. 

It is proposed to determine the structures dynamic and static 
fields of this flow. We will study the influence of the Reynolds 
number and shear-thinning characters on the dynamic field 
and the evolution of head losses. 

The development or the use of numerical calculation software 
is essential to the resolution of such problems. The 
discretization of the system of equations is performed using 
the finite volume method that is implemented in digital 
modelling software "Ansys-Fluent". 

II. Mathematical Model 

The general equations governing the laminar flow of 
incompressible fluids are: 

(i). Continuity Equation 


div(j?P) = 0 (1) 

(ii). Movement Equation 


P~ = -grad(p) 4- divfj) (2) 

-4 _| 

M/feprff 1 is the velocity vector, t the stress tensor, p the 
density and p the pressure. 

The general conservation equations governing the isothermal 
flow in pipes are written by adopting the system of cylindrical 
coordinates. 

In addition the following simplifications are adopted: 
•Stationary flow. 

•Incompressible flow (density is constant). 
•Axisymmetric flow in cylindrical pipe. 

•The radial pressure gradient is negligible. 

(iii). Simplified Continuity Equation 
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3 ini) Sirv) 
dz + dz 


The consistency K of the fluid is found in terms of Reynolds 
number Re and the power index n for a cylindrical pipe and in 
between two plates [4]: 


(iii). Simplified Movement Equation in z Direction 


3l i 3u 13/ 


3lA 3 .. 3u_ / y1 . 


(iv). Simplified Movement Equation in r Direction 


jRe* = S 1 "" f—)" pi4lTCj?TC 

Vlfi + i/ If 


( 11 ) 
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In addition to the following boundary conditions: 


In this chapter, we have established the mathematical partial 
differential equations governing the laminar flow of 
non-Newtonian pseudoplastic fluid and the software used 
"Ansys-Fluent" is responsible for the discretization and 
numerical solution. 


(v) . Condition of non-sliding at the wall 

u (r = R) = v (r = R) = 0 (6) 

(vi) . At the axis of symmetry 

f( r = 0)=^Cr=0) =0 (7) 

3r 3r v ; 

For the problem of the establishment length, we have a 
constant axial velocity profile: 

— (z = 0 j = utP = 1 (8) 

Where r is the pipe radius, u the axial velocity, v the radial 

velocity, is the apparent viscosity and is the mean 
velocity. 

The Simulation was done by the software “fluent” using the 


2. Dynamic Study of Establishment Length 

In this chapter we will study the behavior in the entrance area, 
a non-Newtonian fluid non-complex thermo-dependent 
flowing in a cylindrical industrial pipe, where the axial 
velocity profile constantly evolving until it takes its final 
shape of establishment for laminar flow. 

(i). Flow State in the Pipe Line 

Consider a pipe fixed to a large-sized tank figure (1). At the 
entrance of the pipe, the velocity distribution is flat and 
uniform. During the flow along the pipe, the fluid particles 
adjacent the walls are braked gradually causing the particles 
adjacent to the axis to accelerate so that the flow is conserved. 


BOLOllfpLliW 



Figure 1: Flow state at pipe entrance 

Theoretically for a laminar flow (n=l & Re<2100) the 
establishment length is expressed by [6]: 


model of Herschel-Bulkley [5] for T <T 0 . 

We note that the equation used is Ostwald de Waele, but to — 0,0575^*? (13) 

solve this problem having our study fluid of this type, we 
choose to adopt as law behavior the equation corresponding 

to Herschel-Bulkley fluids in non-Newtonian critical stress. Where B is the pipe diameter . 

Recall that the model of Herschel-Bulkley is given by: 

(ii). Velocity Profile in the Established Flow 


v 


(9) 


Having t 0 “the yield stress threshold” equal to zero the 
preceding equation becomes the Ostwald de Waele equation: 


The velocity profile for an established flow can be defined 
then by: 



(14) 


t = K x l/T L x y 


( 10 ) 


The theoretical velocity profiles for different index n are 
shown in the figure 2: 
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Figure 2: Velocity profile in a cylindrical pipe 

(ii). Geometry and Calculation Procedure 


for the cylinder diameter D = 1 and pipe length L =50D. 





Scaled Residuals 


(iii). Variation of the Establishment length as a function 
of the flow parameters 

From our study and as we fixed the number of Reynolds and 
the structure index n we can find the velocity profile for 
several positions in the pipe as shown in figure 7: 

Velocity Profile 


Figure 3: Studied pipe geometry 

Meshing of the geometry is done in the software “Ansys” and 
a 300x50 mesh was used as shown in the figure 4: 



Figure 7: Evolution of axial velocity for Re=250 & n=0.6 for 
a cylindrical pipe 

From the figure 7 it can be found the establishment length was 
about 13.33 for Re=250 & n=0.6 

The influence of Reynolds number in the establishment length 
was studied by fixing n on 0.6 and varying Re as the results 
shown in the figure 8: 


Figure 4: Typical used mesh 

The study convergence was set for a value of 10e-6 and the 
“QUICK” interpolation and the “SIMPLE” algorithm [8] was 
used as the system configurations and convergence results 



Re=100 
*Re=250 
Re=500 
‘Re=750 


Figure 8: Evolution of the axial velocity for different number 
of Reynolds 

While by fixing Re on 250 and varying n we can find the 
influence of the structure index on the establishment length as 
shown in figure 9: 



n=0.3 
-n=0.4 
n=0.6 
* n=0.8 
*n=l 


ANSYS FLUENT 14.0 (2d. dp. pbns. I 

Figure 6: Results evolution 


Figure 9: Evolution of the axial velocity for different structure 
numbers 
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It is clear that in both cases the axial velocity at the center 
increases gradually as the fluid advance in the pipe, however, 
this development does not occur in the same way in both cases 
studied. 

(iii). Search for Correlations 

By tracing the evolution of the establishment length as a 
function of Reynolds number it can be found that the length 
increases linearly with Reynolds number. 

The results of the simulation are grouped in the table 1 and 
plotted as shown in figure 10. 

Table 1 


Simulation Results 


Reynold 

s 
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D) 
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2.80 
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9.40 
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13.33 

14.22 

14.69 
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19.13 
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26.62 

28.36 

29.40 
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29.23 

34.26 

40.12 

42.46 
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50 



Figure 10: Variation of the establishment length as a function 
of Re for different values of n 

From the results it is shown that the variation of the 
establishment length is directly proportional to the Reynolds 
number, forming a straight line passing by the origin for every 
index n value which can be correlated in the form: 


= uWJe„ (15) 

From our results we can reach to the correlation for the 
cylindrical pipe: 


-— = — + 6„7l1 X Re r (16) 

D 1OT v T2 / n - v J 

The correlation found in the literature [1] gives the 
establishment length for a Newtonian fluid for n = 1, then our 
correlation is in good agreement with [1] with less than 1.7% 
relative error, which validates the numerical results. This 


correlation allows to provide the business flow length for 
various values of the Reynolds number and the index n. 

A same study was done for the case of two plates and the 
correlations were: 




1 / -Q65 

1ZZ \ fi 


+ 2.5?) 


X R& n 


(17) 


III. Singular Pressure Drop 

In this chapter, we will study the singular head losses 
associated with the flow in the entrance area of a cylindrical 
pipe and in between two plates. 

(i). Influence of rheological parameters on the singular 
head losses 

Influence of fluid’s consistency K: To study the influence of 
the consistency K of the fluid on the pressure drops, it returns 
to study the impact of the Reynolds number on them [4]. 



Z{m) 


Re=25 


- Re=50 


Re=100 


- Re =2 50 


Figure 11: Evolution of partial pressure for different Re 
number (n=0.6) 

Corresponding a zero outlet pressure, as noticed the pressure 
drop decreases with increase of Re which is as Re decreases 
the influence of viscosity increases increasing the pressure 
drop. 

Influence of the structure index n: 
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Figure 12: Evolution of partial pressure for different n index 
(Re=100) 
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The singular pressure drop increases with the increase in 
structure index n of the fluid. 


(.LJD') = 2.273a 5 * 3551 x (i?£F T2 )-°- 2135T1+0 - 225fl 


(23) 


(ii). Equivalent Length 


Pressure Influence this 


Pressure drop Js can be found theoretically [2] to be: 


IV. Conclusion 


Is = 


(18) 


This gives the normalized form of the equivalent length: 


D — 64 


(19) 


Where z is the pres sure drop coefficient . 
(ii). Correlations 

Results reached are summarized below: 


The study of the laminar flow of non-Newtonian fluids at the 
entrance of industrial pipes is a part of a complete overall 
study necessary for various singularities due to change in 
directions (elbow) or diameter (expansion, shrinkage, 
obstacles...) for an effective design of heat exchangers. 

Numerical simulations and the use of Ansys-Fluent have 
shown significant efficacy and brought us logical results of 
the laminar flow of pseudoplastic fluids in the inlet region of 
industrial pipes. 

A theoretical approach was used in order to obtain the 
equation for pressure loss and the establishment length. To do 
so, one-dimensional analysis for the fully developed pressure 
drop and the establishment length in the two cases was done. 

Meshes, was performed to study the influence of grid 
refinement on numerical predictions. In addition, a technique 
was used to quantify the numerical results for the total 
pressure drop (Ap). 



— *—0=0.3 


—■— rt=0.4 


—A—0=0.6 


— 0 = 0.8 


After we found the pressure drop coefficient we had made a 
correlation to find the equivalent pipe length and the 
establishment length. 

This study needs continuity and requires the study of the 
influence of various peculiarities already mentioned on 
dynamic and kinetic aspects using the software calculations 
adopted. 
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Figure 13: Evolution of pressure drop coefficient for different 
Re number and structural index n 

The pressure drop coefficient can be correlated from the 
curves in figure 13, and using equation 19 we can find the 
correlated form for the normalized equivalent length. 

For cylindrical pipe: 

= 24.416, S *™ 1 X 

( LJD) = 0 . 3 S 14 a. 5 ' 33 5T1 x {R )- a305T1+0JB76 (21) 

For the two plane plates: 

et= 54.559ff 5ja5B x (JT^r 0313 "- 01 ™ (22) 
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